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ABSTRACT 

The stability of radio millisecond pulsars as celestial clocks allows for the possibil- 
ity to detect and study the properties of gravitational waves (GWs) when the received 
pulses are timed jointly in a "Pulsar Timing Array" (PTA) experiment. Here, we 
investigate the potential of detecting the gravitational wave from individual binary 
black hole systems using PTAs and calculate the accuracy for determining the GW 
properties. This is done in a consistent analysis, which at the same time accounts for 
the measurement of the pulsar distances via the timing parallax. 

We find that, at low redshift, a PTA is able to detect the nano- Hertz GW from 
super massive black hole binary systems with masses of ^ 10® — 10^" Mq less than 
^ 10^ years before the final merger. Binaries with more than ^ 10'^ — 10'* years before 
merger are effectively monochromatic GW, and those with less than ~ 10^ — lO'' years 
before merger may allow us to detect the evolution of binaries. 

For our findings, we derive an analytical expression to describe the accuracy of a 
pulsar distance measurement via timing parallax. We consider five years of bi-weekly 
observations at a precision of 15 ns for close-by 0.5 — 1 kpc) pulsars. Timing twenty 
pulsars would allow us to detect a GW source with an amplitude larger than 5 x 10""'^''. 
We calculate the corresponding GW and binary orbital parameters and their measure- 
ment precision. The accuracy of measuring the binary orbital inclination angle, the 
sky position, and the GW frequency are calculated as functions of the GW amplitude. 
We note that the "pulsar term", which is commonly regarded as noise, is essential for 
obtaining an accurate measurement for the GW source location. 

We also show that utilizing the information encoded in the GW signal passing the 
Earth also increases the accuracy of pulsar distance measurements. If the gravitational 
wave is strong enough, one can achieve sub-parsec distance measurements for nearby 
pulsars with distance less than ^ 0.5 — 1 kpc. 

Key words: pulsar: general — gravitational waves 



1 INTRODUCTION 



Gravitational waves (GWs), ripples in space-time, perturb our four dimensional space-time background, on which the pulsed 
radio radiation from pulsars propagates. Through such space-time background pe rturbations, the GW leaves its fingerprint in 
the a r rival times of pu l sar signals by introducing an extra correlated component ( Sazhin 19781 : Detweiler 19791 : Bertotti et al. 
19831 : IWahlauisj Il987l : [Backer fc Hellingsl 1 19861 ). By timing multiple pulsars (quasi-)simultaneously in a so-called "Pulsar 
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Timing Array" (P TA) experiment, one can extract the GW signal, even in the presence of various unrelated noise contributions 
in the timing data dRomani fc Tavlqijl 19831: iHelhngs fc Downsllliii : iBlandford et al.ll 1984 [Foster fc Backeilll99ol : Ijenet et af 
20041 . 12OO5I : ISesana fc Vecchiolboiol : iFinn fc Lommenlboioh 



In general, GW sources are classified into two groups, the single sources and a stochastic background component. A 
stochastic GW background is a supe rposition of multiple single GW sources, which are inseparable in frequency space 
I Cutler fc Thornell200ll : Ixhornel [l989l ) due to finite frequency resolution of finite data. The techniques of detecting a sin- 
gle GW source and that of detecting a stochastic background are therefore quite different. The detection of a stochastic 
GW background using pulsar timing has already been addre ssed in various studies IPetweilejllOTgl : iHeUings fc Downsllliii : 

Backer fc Hellingdiigsd : Ijenet et alji2005l : lHobbs et aLlboogj ). Here we concentrate on si ngle GW sources. 

Recently, many authors have drawn attention to the single source detection problem Jjenet et al. 20041: Sesana fc Vecchiol 



2OI0I : iBurt et ahlboiol : IPinn fc Lommenl2010l : IVardlev et aLH201ol : lOeng fc Finnl201ol : ICorbin fc Cornishlboioh. because o7 



. . 1) „ „ „ . . , 

2OO9I : iBurke-SpolaojIioioh . In this paper, we focus on the problem of detecting GWs from super-massive black hole binaries 
(SMBHB), with masses from 10® Mq to 10^" Mq and orbital frequencies in the nano-Hertz range, which makes these SMBHBs 
indeed potential GW sources for PTA projects. We will show how PTAs can be used to measure the parameters of SMBHBs 
and their GWs. 

The plan for the paper is as follows. We begin with a qualitative estimation of the parameter space of the detectable 
SMBHB population for PTAs in § [2l We then calculate the parameter estimation error for the GW parameters in § |3] where 
we describe the timing response to a single GW source § 13.11 and present a careful and correct treatment of the corresponding 
statistical problems § 13.21 Finally, we summarize and discuss our results in § ID 



pos- 



sible nearby {z <^1) coale scing binary systems that produce GW signals detectable with PTAs ( Wen et al.ll2009l : Sesana et al 



2 THE DETECTABLE BLACK HOLE BINARY POPULATION 



In order to use the time-of-arrival (TOA) data from a PTA to detect GWs and infer their parameters, two conditions have to 
be met. First, the GW amplitude is large enough such that its statistic is significant to confirm the detection, and secondly, 
the GW frequency is in a frequency window where the pulsar timing technique is sensitive. In this section, we will estimate 
the parameter space of the GW sources, i.e. the parameter space of super massive black hole binaries (SMBHBs), where they 
are observable to practical PTAs. 

For a PTA, one needs at least five parameters to specify its configuration, i.e. the number of pulsars (A^'psr) in the array, 
the accuracy of the TOA (fn), the total observing time span (T), the duration (AT) between two successive observing sessions, 
and the pulsar distances Dpsr- These parameters describe both the ampl itude and frequen cy characteristics of a PTA. For 
detecting a stochastic GW background, we do not need pulsar distances ( Jenet et al. 20051 ). But as we will expl ain later, for 



a single GW source, pulsar distances play an im portant role and need to be taken into account in the analysis IjJenet et al 



2004 



Sesana et ahlboogl : Ipinn fc Lommenlboioh 



T and AT determine the frequency characteristics of PTAs. For example, the frequency resolution of a PTA is A/ ~ l/r, 
and the Nyquist fre quency, the maximal recoverable frequency is given by /max ~ 1/(2 AT ) . In the standard pulsar timing data 



reduction pipeline ( Backer fc Helling^ 19861 : Hobbs et al. 20091 : Lorimer fc Kramer 2005 ) , one fits for the periods and period 
derivatives of pulsars. This removes the lowest-frequency compone nts in the pulsar t iming data, and the minimal recoverable 
frequency /min in the pulsar timing data is /min ~ A/ = l/T (| Jenet et al.l 120051 ). A'^psr and (Jn determine the amplitude 
characteristics of the PTAs. One can detec t a coherent sinusoidal signal among all the pulsars in the PTA, if the amplitude 
of the signal is larger than ( Scargld Il982l ) Cn/-^ NpsrNohs, where A'obs denotes the average number of independent TOA 
measurements for each pulsar. The pulsar distance Dpsi plays two roles in the whole picture. As we will show, it introduces 
another signal component in the pulsar TOA data, i.e. the "pulsar term", w hich increases the signal-to-noise and also gives 
a very long time baseline to investigate the GW evolution (jJenet et al.ll2004l ). 

The mean GW amplitude h from a SMBHB av eraged over the solid angle of all orbital orientations of the binary, in units 
where G = c = 1, is (|Thornelll989l : IWen et ahllioogl ) 



h=^{\h+\2 + \h:,\^) = 2 



\2/3 



\125J 



(1) 



where h+ and hx are the 'plus' and 'cross' gravitational- wave amplitudes, and M, z and Dc are the chirp mass, the cosmological 
redshift, and the co-moving distance of the binary. The ojg is the GW angular frequency at the observer. In the standard 
cosmology model, the co-moving distance is 



dz 



Ho Jo ^!^A + f^M(l + 2)3 

where we used the dark energy density CIa 



(2) 



0.7, matter density SIm = 0.3, and the Hubble constant Hq = 73kms ^ Mpc 
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2.37 X 10 s ^ ( Particle Data Grouj 2008[). The a ngular Irequency relates to the black hole binary chirp mass and the time 



before final merger of the binary via jHugheslboogl ) 



/ 5 s3/8 I 

where im is the time before final coalescence of the binary system at the observer, i.e. the time span between the time of 
'present' to the time of the binary's final merger. The GW radiation takes away the energy of the SMBHB, which induces a 
frequency chirping GW. The angular frequency of the GWs increases at the rate of lj^, where 

'^^-A\256j M5/84V8(i_^^)5/8 • ^' 

The GWs from a black hole binary introduce signals in the pulsar TOA data. It turns out that the GW induced timing sig- 
nal R(t) depends on both the G W strain at the pulsar and the GW strain at the Earth, i.e. R{t) ~ [h{t, 0) — h{t — Dpsr, Dpsr)] /2ajg 



l Estabrook fc Wahlguist 19751 ) q The h{t, 0), the Earth term, is the GW strain at the Earth at the time when we receive the 



pulse; while h{t — Dpsr, Dpsr), the pulsar term, is the GW strain at the pulsar at the time when the pulse was emitted. Thus, 
in principle, the GW from a quasi-circular SMBHB introduces two quasi-monochromatic components in the pulsar timing 
residuals. One of the components comes from the Earth term, the other comes from the pulsar term. Due to the evolution of 
SMBHB system, the frequency of the pulsar term and the Earth term are different. The frequency difference between the two 
components is 5f ~ ujgDpsi/^TT. 

To be observable, the GW induced timing signal has to be large enough compared to the other noises affecting the TOA 
accuracy. For most of the GW single sources, they evolve slowly {5f ^ uug/2TT), thus we can ignore the frequency difference 
between the Earth and the pulsar term and use a single frequency to calculate the amplitude of the induced pulsar timing 
signal, i.e. 

h . a„ 



v/iVpsriVobs 

which is generalized to red noise in Appendix [D] Consequently, 
k5/8 „ n 

A/15/8,1/8 ^ O <7n^c 

We also need to check if the GWs meet the frequency range of PTAs, i.e. fmin ^ cjg/27r ^ /max, which leads to 

StT (1-^^)5/8--'^' - 4^ (1+2)5/8- 

equation. (|6]), the amplitude condition, and the equation. ((7|, the frequency condition, together determine the SMBHB 
parameter space, in which the GWs from such SMBHB will be observable to the PTA. 

It is interesting to see when evolutionary effects become important as discussed by Jenet et al. ( 20041 ) . If the frequency 



difference between the pulsar and the Earth term is larger than the frequency resolution of the PTA, one can detect the 
evolution of the SMBHB. This requires Sf ^ A/ ()SetQi200a '). i.e. 

j.5/8,11/8 ^ 3/5 \3/8 TDpsr , , 

^ ^8^ 1256 J (1 + 2)5/8 • 

It is also interesting to examine when the GWs from SMBHBs are no longer quasi-monochromatic. This is the case, when 
the SMBHB arrives at the dynamical evolution phase, as the orbital radius approaches the innermost stable circular orbit 
(ISCO). The quasi- monochromatic description of the GW is therefore only valid before the binary arrives at the dynamical 
phcise, i.e. the ISCO is a definitive upper bound for the frequency range of the quasi-monochromat ic regime. For an equa l 



mass binary the orbital angular frequency of the ISCO is given by ojorb = /2 ~ 1/(6^''2 2^''5>I) l|Blanchet fc Iveill2003l ). 
Consequently 

tm , 405 

We summarize the results in Figure [T] We can see that PTAs will be sensitive to the SMBHB population with a mass 
range from 10* to 10^" A/q and about 10^ years before the merger. For those SMBHBs with less than ~ 10^ — 10* years, 
the frequency chirp will be visible in the 5-year data. In Figure [1] the non-evolving SMBHBs are believed to be the dominant 
population in terms of numbers, since they have lower masses (10 times smaller) and longer lifetimes (10^ — 10"^ times longer). 
Hence, we will focus on the non-evolving sources in the rest of the paper. The discussion for the evolving GW sources, i.e. the 



^ Here, we ignore the unimportant geometrical factor for this qualitative estimation. 
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Figure 1. The parameter space of SMBHBs as detectable GW sources for a PTA. The right panel shows an enlarged region of the left 
panel. The x-axis is the time span from the 'present' to the final merger in years; the y-axis is the chirp mass of the SMBHB system in 
Solar masses. The corresponding limits are labeled on the line, i.e. 'Amplitude Limit' for equation. JBJ, 'Low Preq. Limit' and 'High Freq. 
Limt' for equation. 10 , and 'ISCO Limit' for equation. ||9ll. The parameter space for evolving and non-evolving GW sources is shaded 
using using light gray and dark gray respectively as indicated in the figure. It should be noted that, although the dark gray region may 
look smaller than the light gray region, the dominant detectable SMBHB population belongs to the dark gray region, where we encounter 
smaller mass and longer rest lifetimes. For our calculations we considered SMBHBs with a and assumed bi-weekly observations of a PTA 
with 20 pulsars, of which the parameters are shown above the panels. 



source in the starred region of Figure [T] will be presented elsewhere. The SMBHBs approach the last stable orbit only about a 
few years before merger, so that for all practical PTA observations, we do not need to consider the dynamical merging phase 
- this is in contrast to the case of ground based GW detectors operating at much higher frequencies. 



3 ABILITY OF A PTA TO MEASURE THE PARAMETERS OF A NON-EVOLVING SMBHB 

In the previous section, we estimated the parameter space for the PTA-detectable SMBHB population. In this section, we 
turn to a detailed analysis of the GW parameter estimation problem. We firstly calculate the pulsar timing response to a 
monochromatic GW, before we determine the statistical error of estimating GW parameters. 



3.1 Pulsar timing signal with a monochromatic GW 



Signals in pulsar TOA data come from many different origins, which include the motion of the Earth, the motion of the 
pulsar (in particular if it is in a binary system), changing relativistic delays in the signal prop agation, the scattering, and 



the dispersing due to space plasma, the gravitational waves in the space time background, etc. l|Manchester fc Taylor 



1977 



Backer fc Helling j Il986l : IStairsI boosl : iLorimer fc Kramerl liooj ) . We are, here, particularly interested in the timing signal 
components related to gravitational waves. 

A persistent GW introduces two different signal components in TOA data, i.e. the "Earth term" and the "pulsar term", 
where the Earth term depends on the GW strain at the Earth and the pulsar term contains the GW strain at the pulsar. 
Due to the phase coherence, the Earth term and the pulsar term interfere with each other depending on the pulsar distance, 
the pulsar direction, and the gravitational wave propagating direction. Because of such pulsar-Earth term interference, the 
pulsar distance also becomes an important variable in the GW detection problem. The most relevant pulsar timing signals 
are then the GW induced signal i?g, the pulsar timing parallax i?par, and the noise n which includes other noise contributions 
to the TOAs. Both GW-induced signal and timing parallax depend on the pulsar distances. The timing residual signal can 
be written as 



R{t) = To + PoiV + iPoPoiV' + Aq 



+ i?g(t) +-Rpar(i) +Jl(t|crn 



(10) 



where the A'' is the sequential pulse number, the To,Po,Po is the initial epoch, period and period derivative for the pulsar, 
while Aq represents the standard Solar system correlations (jLorimer fc Kramen l2005l ). The n{t\a-n) is taken to be white 
zero-mean Gaussian random variables with an RMS level of an to simulate the other noise contributions to the TOA. In our 
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calculation, we fit for Pp , Po , and pulsar sky position using a least-square fitting as in standard pulsar timing pipelines |f| 
i Lorimer fc Kramejiiooil '). 

As shown in Appendix [X] the pulsar timing signal induced by a single GW source is 



R^(t) = feo_ sin(A$/2) 1^^^ cos(2^) + B2 sm{2,p)] cos{ujgt - A'I'/2)(1 + cos^ l) 

^CjkJ ^ _L ' ^ COS \7 

+2 [B2 cos (20) - Bisin(2<j!))]sin(Ljgt- A$/2)cost} , 

where the Bi, B2 and A"l? depend on the geometrical configuration of the pulsar and GW position by 

Bi — (1 + sin /?) cos /3p cos [2(A - Ap)] - sin(2^) sin(2/3p) cos(A - Ap) + (2 - 3 cos^ /3p) cos^ /3 , 
B2 = 2cos/3sin(2/3p)sin(A-Ap)-2sin/3cos^/3pSin[2(A-Ap)], 
A$ = a;gDpsr(l — cos6') . 



(11) 

(12) 
(13) 
(14) 



Here, the A and /3 are the ecliptic longitude and latitude of the GW source, the Ap, /3p are the ecliptic longitude and latitude 
for the pulsar position, the angle t is the orbital inclination of the GW binary source, i.e. the angle between orbital angular 
momentum and the direction to the Earth. The ho = D^'^ {I + zfl'^ is the GW amplitude. 

For a persistent GW source, the pulsar term and the Earth term are coherent and have a phase difference of A$. Such 
pulsar-Earth interference leads to two consequences: 1) it changes the amplitude of the timing signal; 2) it changes the phase 
of the timing signal. One can see the two effects from equation. pi[) . The sin(A$/2) term is the modulation for the signal 
amplitude, while the A"l?/2 in the term lo^t — A$/2 is the signal phase shift (see Appendix 1X1 for the details). 

The pulsar timing response H to the GW signal is the normalized ratio between the amplitude of GW-induced pulsar 
timing signal and the amplitude of the GW, i.e. 

(GW induced amplitude in TOA) x (GW angular frequency) 



H 



where 



lim 

T->oo 



(GW amplitude) 

1/2 



- / Ri{i)dt 



sin(A$/2) 



4(1 



P3 = (P?+P|)(l 



c 2 , 4 
6 cos I + cos ( 



[(Bi - bI) cos(40) + 2B1P2 sin(4(^)] 



(15) 



(16) 



The response 7^ is a function of the GW source position, the pulsar location, and the orientation of the SMBHB orbit. Due 
to the interference between the pulsar term and the Earth term, the response H also depends on the distance of the pulsar 
through the term sin(A$/2) in equation. 1)15^ . Figure [2] shows the timing response as a function of 0, the angle between pulsar 
direction and GW source direction. The response H{Q) is simply the response pattern for a single-pulsar GW detector. As 
shown in Figure [21 the response pattern is made up by a large number of spiky lobes. The envelope of these lobes depends on 
the GW source position and polarization, and the spiky lobes are due to the interference between the pulsar and the Earth 
term. Given the angular frequency of the GW cjg and the pulsar-Earth distance Dpar, one can show that the angular width 
of the each lobe is approximately Siobc ~ l/(2Z)psri^g) ignoring the weak dependence on Q. For practical PTA purposes, the 
i'psr ~ 1 kpc and bj^ ~yr~^, so Siobe ~ 10""^ rad. It is the existence of these spiky lobes which gives PTA experiments the 
power to accurately locate single GW sources, given that the pulsar distances are well measured. 

By measuring the phase difference A$ between the pulsar term and the Earth term and by comparing data from different 
pulsars, the pulsar distance can be inferred. However, the result of such a calculation would have multiple solutions, i.e. we 
cannot discriminate between the case of A$ and the cases of any Atl> + 4n7r, where n = 0, 1, 2, . . ., since we have shown in 
equation. (|11[) that the GW-induced timing signal i?g is a periodical function of Ac[> with a periodicity of 47r. In this way, 
one can only measure the modulo of the pulsar distance by 47rajg'^(l — cosS)"^. On the other hand, if there is another, 
independent way of measuring the pulsar distance, even with a lower but sufficient accuracy to remove such confusion about 
multiple solutions, one can use the GW signal to increase the accuracy of pulsar distance measurements by utilizing the GW 

distance modulo. 

"llellingj 



For millisecond pulsars , such an independent distance measurements is available from the timing parallax ([Backer 



19861: iRvba fc Tavloi 



1991 



Lorimer fc Kramer 



2005 



Verbiest et aklboid '). It is expected that with future radio telescopes. 



such as the SKA ( Kramer fc Stappersll201o[ ] ISmits et al. timing parallax measurements will yield very precise distances 



for nearby pulsars. Indeed, SKA timing parallax measurements promise to be of sufficient accuracy to help us to identify the 
correct solution for A$ and, hence, to measure the pulsar distance with even better precision. 

In order to consider this further, we study a simplified version of the timing parallax term. In reality, the timing parallax 



If the pulsar is in a binary system, the fitting of the orbital motion will also remove some signal power. 
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Figure 2. Response pattern H of a single-pulsar timing response to a single monochromatic GW source. For illustration purposes, the 
pulsar distance is chosen to be small with a value of 200 pc and the GW period is chosen as 5 years, in order to show the structure of 
the response pattern. The GW source is in the 0° angle position, the orbital plane inclination is 90° and the orbital plane coincides with 
the paper plane. In this way the plotted response pattern is, in fact, the term (1 + cos 6) sin(A<I>/2) in equation. | |A7|1 . 



is derived from the detailed motion of the Earth from Solar System dynamics ( Seidelmann 2005h . For the purpose of this 
paper, it is sufficient to keep the leading term of the timing parallax, i.e. assuming a circular motion of the Earth, 



-Rpar(t) — 



3[2(Apsr - A0(t))]cos2/3p 



(17) 



where the term r0 is the average distance between the Sun and the Earth, and A0(t) — 27r(t/lyear) is the ecliptic longitude 
of the Earth at time t. This form of timing parallax assumes that the eccentricity of the Earth orbit is zero. This assumption 
is valid for cases where the pulsar is not too close to the ecliptic poles, i.e. (—89° ^ /? +89°), such that a timing parallax 
signal is not dominated by the Earth orbit's eccentricity. As this will generally be the case, the error of the measured pulsar 
timing parallax distance is (see Appendix [B] for details) 

2 



V Nobs r'i, cos2 



2.34 

3S2 /3ps 



/iVob.y 

V 100 / 



1 kpc 



10 ns 



pc 



(18) 



where A'obs is the number of TOAs. The numerical factor is derived assuming that the time span of pulsar data is longer 
than one year. In a real data analysis, one always uses the full Solar System ephemeris. We compared eq uation. (llSp with 
results from numerical simulations based on TEMPcH and the planetary ephemeris DE405 ( StandishI 19981 ). For pulsars with 
—89° ^ /3 ^ +89°, we find that the simplified version of the timing parallax shown above agrees with the correct result 
derived from TEMPO within a few percent difference, justifying the usage of equation. (|18p for the purpose of the present 
paper. We note that the validity of equation. (|17p comes from the fact that the Earth orbital eccentricity is small and that 
we are investigating measurement accuracies, where the effect of orbital eccentricity is of even higher order. According to 
equation. (118(1 . with a timing accuracy at the 10 to 30-ns level, one can use the timing parallax to measure the pulsar distance 
accurately to a few hght years for pulsar distances of less than 1 kpc. This distance accuracy become comparable to the 
wavelength of the GW, and the timing parallax measurement is therefore indeed a potential technique to remove the pulsar 
distance confusion. Both GW parameters and pulsar distances should thus be estimated from pulsar timing data at the same 
time. In the following, we estimate the corresponding accuracy of the GW parameters and pulsar distances measurements 
based on the signal timing of equation. (|10p . 



3.2 Vector Ziv-Zakai bound for signals with additive white Gaussian noise 

We are, now, going to determine the statistical error of estimating GW parameters using data from a PTA. A well known 
and popular statistical technique to calculate such lower bounds of the statistical accuracies of parameter estimators is the 



^ See |http: / /www.atnf .csiro.au / research / pulsar / tempo / [ 
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Cramer-Rao bound (CRB) based on the Fisher information jFisjllQG^ ). However CRB is also known for p redicting; a too smal l 
value, which can not be achieved in practical cases, especially, when the si Rnal-to-noise ratio (SNR) is low IjChazan et al.lll975l : 
Bell et aljfliii . 1 19961 : [Nicholson fc Vecchiolfliii : IVan Trees fc BellllioOTl '). Bes ides the lo w-SNR problem, the CRB is derived 
in the local sense, which assumes that an unbiased estim ator (UBE) exists ( Fisj 19631 ). Due to this assumption, the CRB 
is not applicable to cases, where the UBE does not exist ( JavnesI 2003h . For example, CRB predicts smaller-than-achievable 
bounds, when there are multiple isolated regions with similar likelihoods in the parameter space. In the present GW parameter 
estimation problem, the pulsar distance confusion introduces multiple equal-likelihood structures for both pulsar distances 
and GW source locations at low SNR. This urges us to use, instead of CRB, other technique s to overcome these difficulties. 
Many statistical bounds were investigated in the p assed decades ( see Van Trees fc Belli ([20071) for a review) . It turns out that 
a group of statistical bounds, which originates from Chazan et al. 1 19751 ). gives more reliable bounds than the CRB; and these 
bounds approach the CRB when the SNR is high and multiple equal-likelihood confusions are resolved. These bounds are 
now referred to as Ziv-Zakai (ZZ) bounds. In this paper, we adopt the ZZ bounds to investigate the error performance of the 
GW parameter estimators. We also calculate the same results using CRB for comparing purposes. 

For a fc-dimensional vector parameter A; (i = l...k), and its estimator A^, the estimation error is quantified by the 
correlation matrix dj = {(A^ — Ai)(Aj — Aj)) . The error for the estimator A^ is just the diagonal part of the correlation matrix, 
Cii. The ZZ bound states that the following inequality holds 



i.e. al 



V < max 



[■Pao,.(A>) + Pao,,(A, + (5A0]P.,min(A„ A, 5A0 dA, 



■qdri , 



(19) 



where Ki can be any constant vector. Following the Einstein summation convention, we sum over the index if it appears 
twice in a single term. The Paq ^(Ai), Paq ii^i + SXi) are the prior probabilities for parameters Ai and Ai -|- S\i given that 
the true value of the parameters are Ao,i. The SXi is a fc-dimensional vector, which maximizes the integral in the rectangular 
bracket under the constrain that HiSXi = 77 in equation. (I19|l . where rj is the scalar parameterizing the integral path. The 
Pe,min(Ai, Ai -|- SXi) is the minimal probability of making an error, when discriminating between the parameter set Ai and 
Xi + 5Xi. The scalar function V is the 'valley-filling function' defined as 



V[/(r,)] =max/(r,-f 7)') 



(20) 



The integration domain is determined by prior information. For example, we may use the fact that the PTA pulsars 
are all located in th e Galaxy, which allows us to confine the integrations for the pulsar distances to the radial size of the 
Galaxy. We refer to iBell et all |l994h for details of the formalism in the above. If Ki is taken to be the fc-dimensional unit 
vector, equation. (|19|l gives the bounds for the parameter estimation errors. We present more details on how we calculate the 
Ziv-Zakai bound for our problem in Appendix [C] 

In order to demonstrate why the ZZ bound is more trustful than classical CRB, we take the GW source position as 
an example. In Figure O we show the Pe,min as a function of GW source location. For the case when the GW amplitude 
is low, there is a high probability of making a mistake in discriminating between the true GW source location and other 
locations in multiple isolated regions. The CRB is flawed by the large local derivative near the most probable value and it 
thus underestimates the error by ignoring other isolated high likelihood regions. The ZZ bound, unlike CRB, is an average 
of the parameter distance 77 weighted by the error probability Pf.min over the entire parameter space. Thus, we only consider 
the ZZ bound to be suitable for a correct treatment of cases with potential multiple solution confusion. 

Since we are interested in calculating the errors of the estimated GW source parameters as functions of the GW and pulsar 
parameters, the priors P^o ; (Xi) and Pxg ^ (Xi+SXi) take forms of Dirac 5-functions, i.e. Pao..(A«) = Pxo_,{Xi+SX^) = <5(A,-Ao,«), 
where Ao,i is the true parameter. With the pulsar signal model specified by equation. (|10[) . one can calculate the ZZ bound for 
the parameter estimators, i.e. for the pulsar distance Dpsr and the GW parameters ^0, '-, A, /3, 0, Ug, (^0. For a PTA with A'^psr 
pulsars, the total number of parameters is A^psr + 7, where the seven represents the GW parameters. Details describing the 
calculation of the ZZ bounds are given in Appendix[Cl We present the results in Figures [H where the errors of the parameter 
estimators are given as functions of the GW amplitude for various configuration of PTAs. For high-SNR cases, the two bounds 
(ZZ and CR) agree with each other within the computational accuracy. For intermediate-SNR cases, the ZZ-bound indicates 
a larger uncertainty than for the CRB predictions. This confirms that the CRB underestimates the parameter errors, which 
is mainly due to the multiple solution confusions as we discussed before. At the low-SNR limit, the ZZ bound converges to 
a fixed value because of prior information about Q,n, i.e. we made use of the fact that the error of the inclination angle t 
and source position error will not be larger than 27r, and that the error of the GW frequency will not be greater than the 
sampling frequency. From the results, one can also see that the CRB underestimates the error of the GW source location and 
GW frequency by one or two orders of magnitudes for intermediate-SNR cases. We conclude that the CRB is ind eed not an 



appro priate tool to investigate the error performance of PTAs as a GW detector, in accordance with the results bv IBell et al 
1 19961 ). 

We present the error of pulsar distance measurements as a function of GW amplitude in Figure [5| As a pulsar distance 
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Figure 3. The minimal probability of incorrectly determining the source location. Here, the GW source is located at the center of the 
white circle. The gray scale shows the minimal probability of making a mistake (Pf.min) when deciding between the true GW source 
location and any other location in the X — /3 parameter space. When the GW amplitude is low (top panel), Pf^min is large (~ 0.5) over 
the whole sky. When the GW amplitude is strong enough (bottom panel), the probability of determining an incorrect GW location is 
lower and the multiple-solution confusion is removed. Thus only the real GW position is left. For intermediate GW amplitudes (middle 
panel), there are multiple isolated positions, which all have high probability of being indistinguishable from the true GW source position. 
For the illustration purpose here, these figures are calculated for a PTA with 40 pulsars, where each pulsar has a 10-ns timing accuracy 
and a distance of 100 pc. The GW frequency is taken to be 1.2 year. The white stars indicate the locations of the PTA pulsars. The 
visible interference pattern is due to the response pattern discussed in the text. 

we use 0.5 — 1 kpc. The measurements of the pulsar timing parallax is for the intermediate-SNR case. Again, the CRB and 
ZZ bounds differ from each other, i.e. CRB underestimates the measurement errors by orders of magnitude. 



4 DISCUSSION AND CONCLUSIONS 

In this paper we have shown that we can use a pulsar timing array to detect the non-evolving single gravitational wave sources 
and measure their physical parameters. One can detect a GW source with a GW amplitude larger than 5 x 10~^^ by timing 
20 pulsars bi-weekly (once every two weeks) for 5 years with a timing accuracy of 15 ns. This agrees with the analytical 
estimation in equation. ([5]). Such observational requirements (i.e. 15 ns timing accuracy) have not been achieved at present, 
but should be possible for future radio-astronomical instruments such as the Square Kilometer Array (SKA) and possibly, for 
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Figure 4. The parameter estimation error as a function of GW amplitude. (Th is t he lower bound for the accuracy in estimating the 
GW amplitude ho, cti is for the orbital inclination angle of the GW source, (ta = ^crj sin"* /9 + tr^ is the lower bound for the accuracy 
in measuring the GW source position, where cr\ and cr^ are the errors for estimating A and /3, respectively. The aui^ is the lower bound 
of the measurement accuracy of the GW frequency. The ho is the GW amplitude. One can see that the error of amplitude estimation is 
not sensitive to the true GW amplitude, while the uncertainties in the other parameters decrease as the GW amplitude becomes larger 
jTakahashi fc Setdliool) . The thick lines are the results from the ZZ bound, and the thin lines are the bounds calculated from the CRB. 
The dashed line is for a GW with a period of 2.2 years, while the solid line is for a GW with a 4.4 year period. The total observation span 
T is 5 years and the duration AT between two successive observation are 2 weeks. The pulsars in the PTA are uniformly distributed 
over the sky at equal distance. Details about the PTA, such as the number of pulsars, timing precision and pulsar distances are given at 
the top of the figure. 



D„ =0.5 kpc D„ =1 kpc D„ =2 kpc 

Psr ^ Psr ^ Psr ^ 




Figure 5. Error for pulsar distance measurements as a function of GW amplitude. With a PTA GW detection, one can measure the 
distance to the used pulsars with improved accuracy. The left, middle, and right panel correspond to pulsars with distances of 0.5 kpc, 
1 kpc, and 2 kpc, respectively. The thick lines (upper two) correspond to the ZZ bound and the thin lines (lower two) correspond to the 
CRBs. The solid line is for GWs with a period of 2.2 years and the dashed lines is for GWs with a period of 4.4 years. It may be counter 
intuitive that the ZZ and CRB agree at the low SNR limit. However, this is the consequence of GW-induced multiple solution confusion 
for pulsar distances as we discussed in the main text. In fact, the difference between ZZ bound and the CRB is due to the GW-induced 
multiple solution confusion of pulsar distance. Because pulsar timing parallax alone has well defined unbiased estimators, the ZZ bound 
and the CRB agree with each other when the amplitude of the GW is small and little confusion is introduced. 
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a limited number of sources, al so for the Five- hundred-meter Aperture Spherical Radio Telescope (FAST) I Nan et al. 20061 : 
Smits et aljboogl : iKramerllioiol ) . For such large telescope sensitivity, radiometer noise will not be the limiting factor to achieve 



a 15 ns-level noise budget. Rather, limiting factors may arise from other physical effects, which we will discuss shortly. 

In our treatment of single GW source detection, we do the novel step of including the pulsar term as an integral part of 
the signal that we are looking for. This is possible since the timing parallax offers the possibility of modelling the pulsar terms 
accurately for nearby sources. We take advantage of the pulsar term in several ways. Firstly, by using the pulsar term, we 
increase the signal-to-noise ratio (SNR). When the pulsar term is completely ignored, the GW signals (only including Earth 
terms) in pulsar timing data are sinusoidal waves with identical phases and the pulsar term becomes a noise source with the 
same amplitude compared to the 'signal' term. The SNR of single pulsar data is then always less than one. Even worse, the 
ignored pulsar-term is sinusoidal signal with unknown phase. It is not a stationary signal. It is only reduced by averaging over 
multiple pulsars but not by time averaging. On the other hand, if we can use the pulsar term as a 'signal' term, we can sum 
the signal from multiple pulsar coherently and the SNR of single pulsar data can be arbitrarily large, depending on the GW 
amplitude. The second important benefit of utilizing the pulsar term is that it allows us to get a precise measurement for the 
GW source location. The interference between the pulsar term and the Earth term introduces spiky lobes in the single-pulsar 
response pattern (see Figure[2)|. In Figure[3l where the small scale ring-type structures come from individual response patterns 
and the large scale quadruple structure is introduced by the coherence of the Earth term. This large scale quadruple structure 
is, in fact, the response pattern, when one only uses the Earth term as a signal. Clearly the small scale spiky structure increases 
the locating ability of PTAs substantially and allows astronomers to search of electromagnetic counterparts of the GW source. 

Thanks to the pulsar timing parallax signal, we are able to include the pu lsar term in our scheme to detect non- 
ev olutionary GW sou r ces. O ur method, first presented in Kramer fc Stappers ( 2O10l ). was complemented recently by the work 
of Corbin fc CornishI (|201GI ) who also include the pulsar term in their analysis pipeline. Their approach to use the pulsar 
term however differs, as they study the chirped signal and ignore the timing parallax signal, while we use the pulsar term in 
combination with the timing parallax signal. 

The pulsar term leads to many potential important applications. For example, Deng fc Finn ( 201GI ) recently point out 
that the curvature of the wavefronts introduces a GW parallax effect, which can be used to measure the SMBHB distance. 
Other possible applications are very interesting to study also. For instance, while the pulsar distance measurement helps the 
GW detection, the GW signal, on the other hand, also helps the pulsar distance measurement. Due to the interference between 
the pulsar term and the Earth term, the GW signal is sensitive to the pulsar distance. In this way, the detected GW signal can 
be used to increase the accuracy of our distance measurements. As we expected, the overall measurement accuracy for pulsar 
distances increases with GW amplitude, as shown in Figure [5l The improvement in accuracy increases faster with respect to 
the GW amplitude for pulsars with smaller distances (Dpsr = 0.5 kpc) than the ones with larger distances (Dpsr = 2kpc). 
This is because the distance measurement using the timing parallax has higher precision for near-by pulsars than for far-away 
pulsars. A smaller number of multiple distance solutions occurs for the nearby pulsars, for which the accuracy of distance 
measurements increases faster when the GW amplitude becomes larger. 

In this paper, we have treated the noise as an additive uncorrelated noise (white noise). Because noise arises from various 
different sources: 1) the instrumental noise, e.g. clock errors, ephemeris errors, polarimetry calibration errors, radio-frequency 
interference, etc.; 2) pulsar intrinsic noise, e.g. micro glitch at even smaller level, switch of spin down rate, pulsar profile 
change, etc.; 3) propagation effects, e.g. DM variation and scattering. Thus the noise is not necessarily white. For a stochastic 
GW background detection using pulsar timing arrays, a co rrect treatment o f the red noise is a very important issue, since 
red noise may destroy any measurable correlation function ( Jenet et al. 2005I '). However the red noise issue is less critical for 
single source detection. This is a consequence of the spectral analysis. Since the signal spectrum for non-evolving SMBH will 
be a single peak at the GW frequency, the SNR of the detectio n is basically the ratio between the GW induced signal power 
and the noise power in the frequency bin ( Yardlev et al. 2O10l ). Thus, red noise may be dominant at lower frequencies, but 
less effective in the frequency bin of the signal from a single GW source. This is fundamentally different from the stochastic 
GW background detection. 

Although several pulsars have timing residuals that are essentially white ( Verbiest et al. 20091 ) at the current level of 
accuracy, it is still an unanswered question as to whether the residuals remain whi te for much higher timing precision or for 
much longer observing span. Se veral authors have modeled the statistics of red noise jFoster fc Backeijll99d :l ICordes fc ShannonI 



2O10l : IShannon fc Cordedboiol ) while some works show the po t ential of improv ing timing accuracy by introducing corrections 



to interstellar medium effects (jllemberger fc Stinebring 20081 : Demorest 20ld ). Indeed, for a few selected pulsars a timing 
precision at a 30~50 ns- level jDemorestlboiol ') is already possible at present, but further studies are needed to see how and if 
this can be improved. 

Beside pulsar timing parallax techniques, one can also determine pulsar distances from other methods, su ch as the Very 
Long Baseline Interferometry (VLBI), pulsar orbital parallax, and relative acceleration ( Smits et al. 2O10l ). Since in this 
paper we mainly focus on extracting GW information and pulsar distances from pulsar timing array data, investigation on 
incorporating these extra information into PTA GW data analysis will be presented in future studies. 

Finally, we have demonstrated how the errors in the parameter estimation can be determined reliably. We have calculated 



Gravitational wave astronomy of single sources with a pulsar timing array 11 



the parameter estimation error using both the CRB and the ZZ bound. Both of the methods give nearly identical results in 
the high SNR region, as we expected. But the two methods deviate from each other when the SNR is low. The CRB is known 
to predict unreachable accuracy in the low-SNR regime, while the ZZ bound gives more trustful results in that region. For 
most applications involving GW detection, the noise contributions dominate above the GW signals, so that the ZZ bound 
turns out to be superior in reliability compared to the CRB. 
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APPENDIX A: PULSAR TIMING RESIDUALS INDUCED BY A MONOCHROMATIC GW 



In this section we consider a non-evolving binary system in a circular orbit as the source of a monochromatic gravitational 
wave. If we choose the coordinate system such that the gravitational wave is i n the g-axis and the ascending node along the 
a;-axis (reference frame of the source), the two polarizations can be written as Jsianchetlbooi l 



h+{t) 



ho cos(ci;gt + $o' 



1 + cos'' t 



hxif) = fto sin(Lijgt + $o) cos t 



(Al) 
(A2) 



<l>o is a constant phase shift and l denotes the angle between the orbital angular momentum of the source and the direction 
from the source to the Earth. The amplitude ho depends on the binary parameters and distance (see Section FS.ll) . For a pulsar 
at distance Dpar, with azimuthal angle tl) and with polar angle 9, the residuals caused by the gravitational wave are given by 
( Estabrook fc Wahlauis3ll975l : lDetweileijll979l : IWahlauistl[l987l ) 



R{t) = (1 + cos 6) [{Hf - Hi) cos(2V') + [H^ - ) sin(27/.)] 
where 



Hi 



Hi 



HI = 



H" 



^sin(A$/2)(l + cos'' l) 



ho 



[cos{ojgt + $0 - A<I>/2) - cos(<I>o - A$/2)] 
sin(A$/2)(2cost) 

[sin(Ljgf + $0 - A$/2) - sin($o - A$/2)] 



(A3) 



(A4) 



(A5) 



where 



A$ = ujgDp^ril - cos 6*) . (A6) 
After dropping constant terms, the timing residuals caused by the gravitational wave reads 

■RW = 7^(1 + cos6i)sin(A$/2) 

[(1 + cos^ i) cos(2i/)) cos{uig^t + <E>o) + 2costsin(2?/>) sin(ajgwt + ^d)] , (A7) 

where "!>' = $0 — A"I>/2. One can easily convert these results to the representation in ecliptic coordinates as used in the main 
text. The representation here has the merit of clearly presenting the physical pictures. For example, it looks like a singularity 
at S = in equation. (|A13|) . but the representation here, equation. (|A7|) . clearly shows that the singularity is canceled by 
geometrical factors of GW polarization. The amplitude modulation (sin(A"I>/2)) and phase modulation (A$/2) can be easily 



understood using the identity e' 



2 sin(A"I>/2)e''"'^^*''^^''''^\ where the interference between two sinusoidal 



signals with a phase difference of A3> introduces an amplitude modulation of sin(A<I>/2) and a phase shift of A$/2 + 7r/2. 

One can also calculate the results using the observer frame. The waveform of GWs from a SMBHB with a circular orbit 
at cosmological redshift z and chirp mass A4 is given by jWahlguistlllisTl : Isiandietlbooel : [llughejboogl ') 



hi 



(A8) 



where the A + and Ax are th e waveforms of the two polarization modes '+' and 'x', respectively. The polarization tensors 



and ef. are (|Wahlauistll 19871 ) 



cos A sin /3 cos /3 
sin A sin /3 cos /3 
- cos^ /3 



sin^ A — cos'^ A sin^ /3 — sin A cos A (sin^ /? + l) 

— sin A cos A (sin^ /3 + l) cos^ A — sin^ A sin^ /3 
cos A sin P cos /3 sin A sin /3 cos /3 

sin(2A) sin /? — cos(2A) sin /3 —sin A cos /? 

— cos(2A) sin /3 — sin(2A) sin /3 cos A cos /3 
— sin A cos P cos A cos /? 

where A and /? are the ecliptic longitude and latitude of the G W source. For a non-evolving SMBHB with circular orbit, the 
GW waveform Ax and A+ of the two polarization modes are (|Wahlauistlll987l ) 

1 , 



(A9) 



(AlO) 



A+ = ho cos t sin(2(;/>) sin(ci;gt) — — (1 + cos t) cos(20) cos(a;gt) 
Ax = —ho cos t cos(2(^) sin(Lijgf) -|- -(1 + cos^ t) sin(2(;/)) cos(a;gt) 



(All) 
(A12) 



The angle t is the orbital inclination of the GW binary source, i.e. the angle between orbital angular momentum and the 
gravitational wave vector. The angle (f> defines the direction of the binary ascending node on the sky. The GW source chirp 
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mass M and the source distance Dc degenerate to the GW amphtude ^o, if one only has the PTA data. Instead of discussing 
the detection statistics using full GW source parameters {M, Dc,uig), it is much better to use the reduced parameter sets 
{ho,ujg), which will be used from now on. 



A GW affects the pu lsar timing by introducing shifts to the observed rotational frequency of the pulsar, u (see refs. lDetweilei 
( 1979h : lLee et al.l l|2008l ). and Appendix [A] for details) 



5v_ 

V 



niUj 



■[h,j{t,0) - h,j(t- 



2 1-COS6' -t^-'-H-zj > (A13) 

where n is the unit vector pointing from the observer to the pulsar and 9 is the angle between the GW source direction 
and the pulsar direction for the observer. The hij{t,0) is the GW strain at the Earth (the introduced Earth term) and the 
hij (t — -Dpsr, Dpsr) is the GW strain at the pulsar (the introduced pulsar term). The explicit forms for n and the cos 9 are 



cos Ap cos /3p 
sin Ap cos /3p 
sin/3p 

= cos /3 cos/?p cos(A — Ap) + sin /3 sin/3p 



(A14) 
(A15) 



where Ap, /3p are again the ecliptic longitude and latitude for the pulsar position. 

The GW indu ced TOA signal iig {t) is the integration of the frequency shift over the observer's local time, _Rg (t) — J dt 
( Jenet et al.ll2005l l. Inserting equation. (jAlljl and equation. (jA12[l into Equation. (|A8[l and integrating equation. (|A13[I over 
time yields Equation, (lllll . 



APPENDIX B: ANALYTICAL CALCULATION FOR THE ERRORS OF TIMING PARALLAX 
MEASUREMENTS 

In this section, we derive an analytic expression for the uncertainty in the distance measurement for pulsars using the timing 
parallax. As a first-order approximation, we assume that the Earth follows a circular orbit. In practice, the Earth orbit is 
not only eccentric, but also shows additional, quite complex deviations from an elliptical orbit due to perturbations by other 
masses in the Solar System. However, because the eccentricity of the Earth orbit is small, and we only need a first order 
estimation of the error in the parallax measurement, our assumption is justified. A more complete discussion on pulsar timin g 
parallax measurements can be found in the literature |Backer fc HeUingj[l98^ : lRvba fc Tavlorlll99ll : lLorimer fc Kramejboosl l. 

We assume that the pulsar TOA signal is composed of the timing parallax term and a white Gaussian noise term 
representing noise contributions from other processes, i.e. 



Rit) = i?par(t) +n(t) 



(Bl) 



Regarding the time series of signals as a vectors, e.g. R — R(ti), i — 1 . . . Nohs, where the A'obs is the number of data points 
and ti is the time of the i-th observation. The probability distribution, /(R|Z)psr), for the TOA data R, with given pulsar 



distance Dpsr, is 



/(R|i3psr) = 



27rcrr, 



Exp 



E 



[R{ti) — Rpariti 

2^1 



Using the CRB, the uncertainty of pulsar distance measurements is 



1 



dln/(R|Dp„) 



(B2) 



(B3) 



/(R|Dpsr)dR 



which corresponds to equation. (|18|) of the main text. We have checked this result against simulations with the pulsar timing 
software TEMPO, which uses the planetary ephemeris DE405 for a highly accurate representation of the Earth orbit in the 
Solar System barycenter. For pulsars which are not too close to the poles of the ecliptic (more than one degree angular 
distance), our analytic equation (|B2[) agrees to better than a few percent. 



APPENDIX C: ZIV-ZAKAI BOUNDS FOR ARRAY SIGNAL WITH ADDITIVE NOISE 

The ZZ bound was first introduced for a scalar par ameter estimatio n problem by Chazan et al. ( 19751 ). a nd was later extended 
to gen eral vector parameter estimation problems ( Bell et al. 1994h . Its first use for GW studies was by Nicholson fc Vecchiol 
( 1998^ ). [ n this paper, we only summarize the key results of the ZZ bound and refer readers to the literatures cited for further 
information. 
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We suppose that we are deahng with timing data from an array of A'^psr pulsars. The data of each pulsar is a A^'-point 
TOA serial signal. We denote the whole timing data set by Rj , where the superscript is the pulsar index, and the subscript 
is the data index. For example, the R^j is the j-th data point of the i-th pulsar TOA. The timing data is just a A'psr-times-A'^ 
dimensional vector R and their components R) takes the form of equation. (|10|l . The signal part of R is completely determined 
by the fc-dimensional parameter Xk (GW parameters plus pulsar parameters). The noise part is a random variable following 
independent Gaussian statistics. We have 



R'j = s}(Afc) + n), 



(CI) 

2 {i) 



where the signal part s = Rg+Rpar, and the noise n is assumed to be an un-correlated white Gaussian noise, i.e. (wj?!*) = a 
and {n^jin^j ) = for any j' 7^ j or i' 7^ i. The cri'' is the RMS level for the i-th pulsar TOA noise from other contributions. 
The P£,min(Afc, X'k) is the minimal probability of making a mistake in deciding which pa r amete r set, the A^ or X'f., to choose, 
when giving data R. For the case where A and A' have equal prior, as shown by Kassam 1 19881 ), the following likelihood ratio 
test minimizes the probability of making mistakes. 



choose 



\, if /(R-I^fc) 

^fe' " /(R|a;) 

y -f /(R|Afc) 

^fe' /(R|A' ) 



> 1 
< 1 



(C2) 



where the /(RjA^) is the probability distribution function of the data R, when the parameters take value of Afc. Following 
the likelihood test in equation. (IC2|) . one can show 



p!,min(Afc, Aj.) 



/(R|Afc) 
/(R|A'J 



> 1 



/(R|Afc) 

/(R|A'J 



< 1 



Aj 



For the uncorrelated white Gaussian noise, the likelihood ratio /(Rj Afc)//(RjA^.) is 



/(R|Afc) 




Exp 




(^K^o..) 


- si{Xk) + nlY 


/(2aS«)" 






Exp 






-sj(A;) + n5)' 


/(2a^;«)" 



(C3) 



(C4) 



Integrating over the random vectors n* , one derives 



/(R|Afc 



> 1 



Ao,i 



1 



Erfc 



s){X'k)) 



and the P 



< 1 



Ao,fe ) can be evaluated in very similar fashion. The minimal error probability Pt.min is then 



(C5) 



where 



2Erf(C), 



JVpsr 

E 



E 



N 

3 = 1 



S){Xk) 



2 a- 



2 (i) 



(C6) 



(C7) 



equation. (|C7|) and (|C6[) show that the larger difference between the signal of A^ and the signal of AJ,, the larger ( is, and the 
less probability one will make a mistake in discriminating between the two signals. 



APPENDIX D: GENERALIZATION OF THE AMPLITUDE CONDITION TO THE RED NOISE CASE 

In the main text, equation ^ is valid only for white-nosie. Since equation ([5]) is quite useful to estimate the detectable 
amplitude of monochromatic GWs, we generalize it here also for the red-noise case. The easiest way of deriving equation ([5]) 
is by using a spectral analysis in frequency space. 

We focus on the frequency bin, which contains the GW signal, i.e. the frequency band of [/ — §A/, / + | A/]. For signals 
from A'^psr pulsars, the total GW signal power Ps in the frequency bin is 

and the noise power is 

Pn ^ jVp,,g„(/)A.f = ^P-^"(/)/"-- , (D2) 

JVobs 
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where the Sn{f) is the spectra of n(t), the noise component. The /max, defined in Section [2j is the bandwidth of the signal. 
The detectable condition is then Pa > Pn, which leads to 



^ y ■S'n(/)/max fDSl 

2ojs ^ V ^P-iVobs ■ ^ ' 

In this way, the generalization of equation ([5} to the red-noise case is simply to replace the RMS noise level an with the value 
of cTe = \/ Sn(.f)/max, the effective level of noise. 



